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Relations between the thermodynamics and dynamics of supercooled liquids approaching a glass 
transition have been proposed over many years. The potential energy surface of model liquids has 
been increasingly studied since it provides a connection between the configurational component of 
the partition function on one hand, and the system dynamics on the other. This connection is 
most obvious at low temperatures, where the motion of the system can be partitioned into vibra- 
tions within a basin of attraction and infrequent inter-basin transitions. In this work, we present 
a description of the potential energy surface properties of supercooled liquid water. The dynamics 
of this model has been studied in great details in the last years. Specifically, we locate the min- 
ima sampled by the liquid by "quenches" from equilibrium configurations generated via molecular 
dynamics simulations. We calculate the temperature and density dependence of the basin energy, 
degeneracy, and shape. The temperature dependence of the energy of the minima is qualitatively 
similar to simple liquids, but has anomalous density dependence. The unusual density dependence 
is also refiected in the configurational entropy, the thermodynamic measure of degeneracy. Finally, 
we study the structure of simulated water at the minima, which provides insight on the progressive 
tetrahedral ordering of the liquid on cooling. 

PACS numbers: 61.43.Fs, 64.70.Pf, 66.10.Cb 
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I. INTRODUCTION 

In recent years, numerical study of model liquids in 
supercooled states has been helpful to clarify the physics 
of the glass transition The availability of long tra- 
jectories in phase space offers the possibility of closely 
examining the changes in supercooled states that are re- 
sponsible for slowing down the dynamics by 15 decades in 
a narrow temperature range approaching the glass transi- 
tion. Although current computational studies are limited 
to times shorter than f« 1 /xs (as opposed to real liquids, 
for which the dynamics can be studied up to « 10^ s), 
a coherent picture of the glass transition phenomenon is 
beginning to emerge. 

In addition to characterizing changes in the dynam- 
ics, recent studies have demonstrated the utility of ex- 
amining the underlying potential energy surface (PES) 
as an aid to understanding the properties of supercooled 
liquids connecting the dynamics to the thermodynam- 
ics and the topology of configuration space (2|-||]. This 
connection is most obvious at low temperatures, where 
the motion of the system can be partitioned into motion 
confined within a single potential energy basin with infre- 
quent inter-basin transitions. It has been shown that at 
sufficiently high temperature (at constant volume), the 



system explores always the same distribution of basins, 
and that the average basin energy is nearly temperature 
independent. Below a crossover temperature — which is 
coincident with the onset of a two-step relaxation in the 
decay of density fluctuations |^ — the system starts to 
populate basins of progressively lower energy, but which 
are less numerous. 

The thermodynamic approach based on the analysis 
of the PES, following the formalism proposed by Still- 
inger and Weber [l0|] , has become a powerful formalism 
for the interpretation of numerical data, both in equilib- 
rium [^-^ and in out-of-equilibrium conditions JTT| , p^ . 
The degeneracy of the energy minima, i.e., the number 
of basins with a selected minimum energy, has been quan- 
tified for several model systems and used to calculate the 
configurational entropy S'conf , from which an "ideal" glass 
transition (in the sense of Kauzmann, Adam, Gibbs, and 
DiMarzio |3|-^) has been estimated 

In this paper we present a detailed investigation of the 
properties of local potential energy minima, or "inherent 
structures" (IS), sampled by the extended simple point 
change (SPC/E) model of water |l^]. Previous studies on 
the PES for models of water have shown the relevance of 
this approach to the deepening of our understanding of 
both structural and dynamical properties of liquid wa- 
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ter [[7|,|17H23|]. Specifically, we calculate the temperature 
and density dependence of the basin energy over a wide 
range of temperatures and densities. We also study in 
detail the shape of the basins in configuration space and 
estimate their degeneracy. The information presented 
provides a detailed characterization of the PES and fur- 
nishes all information required to explicitly write the liq- 
uid Helmholtz free energy in a wide temperature and 
density range. Finally, we study the geometrical arrange- 
ment of the molecules at IS minima to better understand 
the changes that take place in the liquid on cooling. Fo- 
cusing on the IS allows us to eliminate thermal effects 
which complicate the temperature dependence. In par- 
ticular, we focus on the fraction of water molecules that 
are four-coordinated. 



II. SIMULATIONS 



The majority of the state points studied here are from 
the molecular dynamics simulations of the SPC/E model 
performed in Ref. ||2^. The simulation methods are dis- 
cussed in Ref. |24|. We have performed additional simu- 
lations of ice Ih, so that we can compare the IS properties 
of the crystal with those of the liquid. The simulations 
of ice consist of a periodic box containing 432 molecules 
with dimensions 2.634 nm x 2.281 nm x 2.151 nm for den- 
sity p = 1.0 g/cm^. The dimensions are uniformly scaled 
in order to obtain other densities. The box dimensions 
have been optimized to generate the lowest energy con- 
figuration at density 1.0 g/cm"^. Proton disorder in the 
initial configuration is generated by identifying closed hy- 
drogen bond loops, and exchanging hydrogens between 
molecules, as described in Ref. p5| . 

The ice simulations have been performed at p = 0.90, 
0.95, 1.00, and 1.05 g/cm^ and temperature T = 194 K. 
The thermodynamic properties are summarized in Ta- 
ble |. The equilibration time for these sample is far less 
than that of the liquid at the same temperature since only 
the vibrational degrees of freedom need to be relaxed. 

We have performed conjugate gradient minimizations 
p6t to locate local minima on the PES closest to any 
given instantaneous configuration. We use a tolerance of 
10-15 kJ/mol in the total energy for the minimization. 
For each state point, we quench at least 100 configura- 
tions taken from two independent trajectories. While 
each configuration is not necessarily separated by the 
typical relaxation time of the system, the set of points 
quenched typically spans « 25 times the relaxation time 
of the intermediate scattering function. 



III. INHERENT STRUCTURE PROPERTIES OF 
SPC/E 



A. Thermodynamics in the Inherent Structure 
Formalism 

Stillinger and Weber formalized the concept of a basin 
in the potential energy surface by introducing the in- 
herent structure formalism [|o|. The set of points that 
map to the same minimum via steepest descent are those 
which constitute a basin, and the minimum of a basin 
is the IS. This approach is particularly well suited to 
simulated liquids, since it is possible to explicitly cal- 
culate the steepest descent trajectory from an equilib- 
rium state point. Moreover, the partition function Z 
can be explicitly written in terms of the basins. In the 
isochoric-isothermal (NVT) ensemble, for a system of N 
rigid molecules 



Z = \ 



j cxp{-V/kBT)d'^i 



(1) 



which can be written as a sum over all basins in config- 
urational space, i.e. 



Z = X 



-6N 



^ expi-eis/ksT) x 



(2) 



exp{~{V^eis)/kBT)d^r 



basin 

Here A = h{2TTrakBT)~^/^ is the de Broglie wavelength, 
V = V{y'^) is the potential energy as a function of the 
atomic coordinates, e/5 is the energy of the IS and Rbasin 
is the configuration space associated to a specific basin. 
The model system we consider here, namely SPC/E wa- 
ter, has six degrees of freedom for each molecule. It is 
natural to introduce f2(e/s), the number of minima with 
energy e/s, and the free energy of a basin with basin 
energy e/s /(T, e^s) (the "basin free energy"). 



f{T,eis) = -ksTln 



1 



-6N , 



(3) 



basins ^ba 



eM-{V ^eis)/kBT]d^r\ , 



The asterisk denotes the fact that the sum is constrained 
to basins of energy ejs- Eq. (^) accounts for both the 
basin structure surrounding the minimum and the kinetic 
degrees of freedom. The complete partition function can 
be written as the sum over all possible ejs values. 



Z = ^ n{eis) exp 



e/s + f{T, eis) 



or 



Z = ^ exp 



T5'conf(ejs) + eis + /(T, ejg) 
kBT 



(4) 



(5) 



where the configurational entropy 5'conf(e/s) = 
/cb ln(r2(e/s)). In the thermodynamic limit, the cor- 
responding Helmholtz free energy F{T, V) is given by 
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F{V,T) = EisiT) + f[T,Eis{T)) - T5co„f(i?/s(r)) 

(6) 

where Eis{T) is the thermodynamic average of e/5 and 
solves 

dF[V, T)/deis = 1 + df{T, eis)/deis - (7) 
TdS 

colli 

(T, eis)/deis = 0. 

Eis{T) can be numerically calculated by estimating the 
IS which are populated by a system in equilibrium at 
temperature T and fixed volume V . Hence, if a good 
model for /(T, e/s) is available, then the e/s dependence 
of Sccmi along isochores can be estimated by integrating 
Eq. (13). Note that Eq. (0) shows that, if the basin free 
energy does not depend on e/5, then the configurational 
entropy is the only quantity controlling the T-dependence 
of Eis- In other words, the statistical mechanics of the 
basins completely decouples from the vibrational dynani- 



Sconi{Eis) can also be calculated by studying the prob- 
ability distribution P{Ejg,T), i.e. the probability that 
the liquid — in equilibrium at temperature T — populates 
the inherent structure Ejs- Indeed, from Eq. (pf) 



sho w the Kauzmann energy Ek, 
sec. 



which we discuss in 



HID 



'S'conf(e/s) = ks In P(e/5,T) + ejs/T + 
fieis,T) + kB InZ(r). 



(8) 



Hence, Eq. ^ gives Sconf up to an unknown constant 
kBlnZ{T). This "histogram technique" has been re- 
cently used to estimate the configurational entropy for 
a binary- mixture Lennard- Jones system [^J^,p^ . 



B. Inherent Structure Energy as a function of p and 

T 



The p- and T-dependence of Eis for the SPC/E poten- 
tial for a more limited range of T was recently reported 
in Ref. J22]. Our results for the IS energies are shown 
in Fig. |VI|(a) as a function of p and in Fig. VI (b) as a 
function of T. 

The IS energy of the liquid along isotherms shows hints 
of negative curvature for T < 230 K (Fig. 0(a)). The 



presence of this negative curvature can also be observed 
in the instantaneous equilibrium configurations [ pO| . 
This curvature yields a negative contribution to com- 
pressibility @ if^i = V[{d^U/dV^]T -T{d^S/dV^)T], 
which might be related to a low-temperature critical 
point in SPC/E Fig. VI also shows the IS energy for 
ice Ih (which coincides with the ground state energy). At 
the lowest T studied, Ejs of the disordered liquid is still 
significantly greater than that of ice . Note that ice 
is not the thermodynamically stable crystalline form for 
SPC/E Hi, and the stable SPC/E crystal would most 
likely have a lower ground state energy. However, the 
stable SPC/E crystal does not correspond to any of the 
experimentally-known forms of ice p4|. Fig V](a) also 



Fig. VI (b) shows that at high T, Ejs is nearly T inde- 
pendent. For T < 350 — 400K, there is a rapid decrease 
of Ejs with a weak density-dependence. At low T, Eis 
depends linearly on 1/T, as shown in Fig. 0(c). The 
1/T dependence of Eis can be derived from the parti- 
tion function provided ^{eis) is Gaussian, and that the 
basin free energy does not depend on e/s Ppql. Further- 
more, as found for the Lennard-Jones liquids^, the T at 
which the IS energy starts to decrease correlates with the 
temperature at which a two-step relaxation starts to be 
observed in all characteristic correlation functions (see, 
e.g.. Fig. 13 of Ref. [||) 



C. The Basin Free Energy: Density of States 



We next focus on the shape of the IS basin with the aim 
of developing a model for the basin free energy. In the 
harmonic approximation, the basin free energy is given 
by 



6W-3 



F{Eis,T) = keT Mf^'^^/kBT)] 



(9) 



the free energy of a harmonic oscillator with frequency 
spectrum uii. The values ujf are the eigenvalues of the 
Hessian matrix, defined by the second derivative of the 
potential energy with respect to the molecular degrees of 
freedom at the basin minimum. The mass and moments 
of inertia of a molecule have also been absorbed in the 
definition of LUi. The distribution of uji, called the density 
of states (DOS), is shown in Fig. |v^(a) for three different 
state points at T = 210 K. The pronounced minimum 
at « 400 cm^^ separates the translational modes (at 
lower frequencies) from the rotational modes (at higher 
frequencies). At larger p, the peaks in the DOS broaden 
due to the disruption of the H-bond network, which we 
will discuss in Sec. [V, For comparison, we also show the 
DOS for ice Ih, where we see a clear separation between 
the translational and rotational modes. 

The normal mode spectrum contributes to the basin 
free energy via the term 6{ln(huj)) (per molecule), where 
the brackets denote an average over the DOS and over 
different configurations. The dependencies on both T and 
p of (hi{huj)) arc shown in Fig. ^(b). The dependence 
of 6fcB(ln(?iLj)) on Eis is shown in Fig. VI(c). The av- 
erage basin frequency is larger in deeper basins, showing 
that the basins become increasingly "sharp" on cooling. 
This is in contrast with the Lennard Jones case, where 
the basins become broader on cooling 1 1 1 jlj] . The aver- 



age curvature of the IS basin at high density has a weaker 
i?/s-dependence (and hence T-dependence) than those at 
low density, but are generally larger than the curvature 
at low density. 
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Fig. V] shows the harmonic free energy estimate of 
Eq. as a function of p. The range of values of the vi- 
brational free energy of Fig. VI (6 to 9 kJ/ mol) is not very 
different from the range of values of Eis of Fig. [v|(a) 
(from —55 to —60 kJ/mol); thus both make a significant 
contribution to the free energy of Eq. (||) . 

The basin free energy estimated using the harmonic ap- 
proximation of Eq. (|9[) can be used as a starting point for 
estimating the true basin free energy. For a more precise 
quantification of the basin free energy, we must consider 
anharmonic contributions to the free energy. Indeed the 
basins of the SPC/E are anharmonic. This can easily 
be seen by considering the difference C/vib = U — Ejs', 
for a molecular system in the harmonic approximation, 
[/^ib = {6/2)kBT. Fig.^shows a marked deviation from 
harmonic behavior. However, as we discuss next, the 
anharmonicity does not have a strong Ejs dependence. 
Development of techniques for the estimation of the an- 
harmonic contribution to the basin free energy would be 
very useful. 



D. Basin Degeneracy 



may be available to this system ||2^ , provided the Ek es- 
timates are not significantly affected by unknown errors 
in the extrapolation procedure or in the reference value 
for S'conf(-Bo) used. 

To confirm independently the validity of the approach 
followed to estimate Ejs dependence of Sconi, we com- 
pare the value obtained from Eq. ( p^ ) with the value 
obtained using Eq. (^ in Fig. VI (c). The agreement be- 
tween the curves is remarkable. Moreover, the overlap for 
different P(e/s,T) distributions after the harmonic Ejs 
dependence is taken into account, suggests that there are 
no other systematic .E/s-dependent contributions. Any 
remaining T-dependent contributions are absorbed in the 
unknown Z{T) function. 

The results reported in this section provide a detailed 
analysis of the inherent structures and of their basins. 
This information can be used to develop a detailed free 
energy expression for the SPC/E Calculation of the 
T dependence of ^conf was presented in Ref. ||^ , to probe 
the relation between the configurational contributions to 
thermodynamic quantities and the liquid dynamics in su- 
percooled states. 



A key element in the description of the configuration 
space in the IS thermodynamic formalism is 0(6/5 ), the 
number of basins with energy e/g. The correspond- 
ing configurational entropy Sconi — i-e. the logarithm of 
^{eis) — can be calculated by integrating Eq. (|^). Using 
the harmonic approximation of Eq. (joj) for the basin free 
energy, we obtain 



Sconi{Els) — Sconf{Eo) 



dE 



IS 



T 



kB{Huj{Ejs)/^{{Eo))). 



(10) 
(11) 



Fig. |V|(a) separately shows the contribution from 
J dEjs/T, and the contribution associated with 
(ln(a;(i?/5')/w((£'o))}- The harmonic contribution is not 
negligible. To obtain the Sconf in absolute scale, the value 
of S'conf at the reference point Eq is needed. We have used 
reference values obtained in Ref. [Q by independently 
calculating the absolute value of the entropy via thermo- 
dynamic integration from the known ideal gas reference 
point. 

The complete results for all densities are shown in 
Fig. VI (b). In the same graph we show a fit to S'conf (^^/s) 
using the form 

SconfiEis) - AiEis - EkY + B{Eis - Ek) (12) 

where Ek, the Kauzmann energy, is the Eis value at 
which the configurational entropy vanishes. This func- 
tional form is equivalent to the Gaussian distribution of 
Vl{eis) Q. The resulting Ejs values, which provide an 
indication of the p dependence of Ek, are reported in 
Fig. VI (a) . The Ek values suggest that disordered states 
with energies comparable to the ordered crystalline states 



IV. STRUCTURE 



On lowering T the dramatic changes in the IS ener- 
gies are known to be accompanied by equally dramatic 
changes in the dynamic properties of the instantaneous 
configurations. However, examination of simple measures 
of structure of the instantaneous configurations (such as 
the pair distribution function) do not reveal such obvi- 
ous changes; rather, there is a very gradual change, with 
the structure becoming slowly more well defined. A more 
careful analysis of structure is is required to see signifi- 
In the following, we focus our atten- 



cant changes |35 



tion on the structural changes that can be observed by 
studying the IS to try to obtain a more clear picture of 
the structural evolution of the system on cooling. The 
results we present are complementary to the results re- 
cently reported for the same system along similar lines 
of thought [E2| and expand on previous work pGjBq]. 



A. Pair distribution function 



Fig. ^ shows the oxygen-oxygen pair correlation func- 
tion for both the equilibrated liquid and for the inherent 
structures at various T and p. We first focus our atten- 
tion on the T dependence along the p = 1.00 g/car' iso- 
chore. For equilibrated configurations as well as inherent 
structures, the first and second peaks of the pair correla- 
tion function become better defined upon decreasing the 
temperature. Further, there is a systematic reduction 
of the intensity between the first and second neighbor 
peaks. At higher density (p = 1.40 g/cva?), the behavior 
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is somewhat different; for the instantaneous configura- 
tions, there is a clear T dependence, while for the IS, the 
T dependence is nearly negligible. For comparison, the 
model atomic liquid studied in |^,0[ , shows virtually no T 
dependence of g{r) in the IS (FigM. Hence the behavior 
at large p is more like that expected for a simple liquid, 
consistent with the disappearance at large p of many of 
the anomalies that differentiate water at ambient density. 

We also show the behavior of g{r) for various p at 
T — 210 K in the bottom panels of Fig. VI. On increas- 
ing p, both the instantaneous and inherent structure con- 
figurations show an increase in intensity at r sa 0.32 nm, 
a trait already known to develop due to distortion (and 
eventual interpenetration) of the hydrogen bond network. 
For water, it is known that the preferred nearest-neighbor 
geometry is tetrahedral. This tetrahedral ordering is 
obvious at low densities, where peaks at 0.28 nm and 
0.45 nm are the expected distances for a perfect tetra- 
hedral lattice with first neighbor separation of 0.28 nm. 
To see more clearly the tetrahedral nature of the liquid 
at all these temperatures and densities, we focus on the 
neighbor statistics at the various state points. 



B. Neighbor Changes 

We first consider the the average number of neighbors 
that a molecule has within a sphere of radius r, which is 
obtained from integration of g{r): 



n{r) = Anp / r'^ g{r')dr' . 



(13) 



We show n(r) in Fig. ^ as a function of p at the lowest T 
studied. We see a plateau in nir) almost exactly equal 
to four for all densities. Therefore, even at large p, the 
liquid has short range tetrahedral order. However, the 
rapid growth of n{r) at large p highlights the distortion 
and interpenetration. 

To quantify the T dependence of the structural changes 
we calculate the distribution of the number of neigh- 
bors a molecule has within a distance 0.31 nm (arrow 
in Fig. H), roughly corresponding to the first minimum 
in g{r) at low density, shown in Fig. VI. We choose the 



first minimum in g{r) at low density to emphasize the 
tetrahedrality of the liquid; for ice, all molecules would 
have 4 neighbors, while in high density liquid configu- 
rations, the distortion leads to a significant number of 
molecules having more or less than 4 neighbors in the 
first shell. Fig. || shows the histograms of the fraction of 
molecules with a given coordination number. At low den- 
sity, the histogram for the instantaneous configurations 
changes from a rather broad one to one that is peaked 
around the value 4 as the temperature decreases, as we 
expect. The same trend is visible for the inherent struc- 
tures, although even at high temperatures, the distribu- 
tion is quite narrowly peaked around the value 4. Such 
a comparison permits us to make a separation between 



deviations from four-coordination arising from thermal 
agitation, and that arising from configurational change. 
This less marked T dependence of the IS can be seen for 
all p. At larger p, the distribution is still peaked around 
4, but there is significant fraction of molecules with are 
not four-bonded — these are molecules are most likely 
participating in the so called bifurcated bonds pO| . 

To quantify the changes in the tetrahedrality as a 
function of T for each p, we plot the fraction of four- 
bonded molecules fi as a function of T in Fig. |l^. Those 
molecules which are not four bonded represent the set 
of bifurcated bonds j2^. At low density, it is interest- 
ing to note that for both the equilibrated configurations 
and inherent structures, this /4 is close to 1 at the low- 
est T simulated. Indeed, for T somewhat lower that 
we can currently simulate, it appears that all molecules 
would be four bonded, and hence form a perfect ran- 
dom tetrahedral network. A simple extrapolation of /4 
for the lowest densities, displayed in Fig. |ll|, shows that 
r(/4 — 1) appears slightly lower than the mode-coupling 
transition temperature Tmct [ p4|j3^j3^ but well above 
Tk- Additionally, it appears that both the instantaneous 
and IS configurations appear to reach a random tetra- 
hedral network at roughly the same temperature. The 
close correspondence of T{fi = 1) and Tmct suggests 
that this crossover in structural change may be the con- 
trolling mechanism for the crossover in dynamic proper- 
ties at Tmct [p| j39| , |40| . However, the extrapolation does 
not allow us to unambiguously associate these temper- 
atures. Since the tetrahedral geometry is the "ideal" 
configuration at these low densities, there are unlikely 
to be any significant structural changes in the liquid for 
T < T{fi = 1); to the extent that there is further struc- 
tural change for T < T{f4 = 1), the rate of change must 
be significantly different than for T > T(/4 = 1). Anal- 
ysis of the neighbor statistics in the IS at lower T might 
shed some light on the hypothesized change in the dy- 
namics from that of a fragile liquid to a strong liquid 
[|l|-|3|, a widely debated topic Q. 

At larger densities increases, it is apparent that a ran- 
dom tetrahedral network is never reached. This is ap- 
parent in the fact that there are far fewer four-bonded 
molecules in the first neighbor shell at the higher densi- 
ties. Hence, no sharp change in the structural develop- 
ment as a function of T is expected. Consequently any 
crossover from fragile to strong behavior must become 
less pronounced. 



V. CONCLUSIONS 



We have characterized the properties of the PES basins 
in the configuration space of the SPC/E model over a 
range of densities in the supercooled regime. We have 
shown that the T dependence of Ejs is qualitatively very 
similar to the previous observations of simple liquids. We 
have also shown the importance of accounting for the 
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shape of the basins when characterizing the thermody- 
namic properties of the IS subsystem. In particular, we 
found the there are significant T dependent anharmonic 
effects. This detailed information on the basins should 
be useful for any future studies that focus on the dynam- 
ics of exploring the PES, or the thermodynamics at low 
T. Such a classification of basin properties would also be 
useful for other model glass formers, and might help to 
highlight the differences between fragile (such as OTP) 
and strong (such as Si02) glass formers [ [l5[ . 

In addition, we presented results for the structural 
changes of the IS; these changes are more pronounced 
than what can be seen from equilibrium configurations. 
In particular, the progression of the structure toward a 
random tetrahedral network at low T holds promise for 
a more physical and intuitive understanding of the glass 
transition in water. Below the temperature of crossover 
to a tetrahedral network, no further structural arrange- 
ment is expected, and the configurational entropy of the 
liquid may be nearly "frozen in" at the value that cor- 
responds to the random tetrahedral network (plus the 
residual defects that would be present at concentrations 
varying with density) . Thus the rate of change in entropy 
of the liquid may be expected change substantially near 
the crossover temperature, resulting in a significantly 
lower value of Tk, compared to the value that may be 
expected if the rate of change above the crossover temper- 
ature persisted. Similarly, because of the significant tem- 
perature dependence of the fraction of bifurcated bonds 
- which facilitate structural rearrangement - above the 
cross-over temperature, and the relative constancy below, 
the temperature dependence of the dynamical properties 
may show a corresponding crossover. 

Finally, we call attention on the possibility of study- 
ing physical aging in this model-system, starting from 
the thermodynamic description of this system. We note 
that, at p — 1.40 g/cm^, there is almost no variation of 
the basin curvature (Fig. [v| ) on Ejs , nor does the struc- 
ture change significantly (Fig. [V^ . Hence these high den- 
sity state points may offer an ideal opportunity to check 
if the out-of-equilibrium dynamics at very low T can be 
still related to equilibrium states of the system . 
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TABLE I. Summary of the thermodynamics properties 
obtained from simulations of Ice Ih, all at T = 194 K. 
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FIG. 1. Inherent structure energy, (a) Density depen- 
dence of the Eis for hquid configurations, and ice Ih. The 
bottom curve is Ek{p), obtained by fitting S^onf using 
Eq. (p^). Note that the three lowest isotherms, as well as 
Ek{p) display inflections. The solid lines are shown as a 
guide to the eye. (b) T dependence of Eis for selected den- 
sities. For T < 400 K, Ejs starts to decrease rapidly, cor- 
responding to the fact that the systems populate basins of 
lower energy, (c) 1/T vs Ejs for all studied densities. The 
linear relationship at low T support the possibility that the 
eis values are Gaussian distributed. Curves for different 
densities have been shifted by multiplies of 0.2 for clarity. 
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FIG. 2. Shape of the basins surrounding the inherent 
structure, (a) The DOS of the liquid at T = 210 K. The inset 
shows the DOS for ice Ih at p = 1.00 g/cm^ for comparison, 
(b) Change of harmonicity of the basins (c) (log(7itj)) vs Ejs 
for different densities. 
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FIG. 4. Anharmonicity of the basins as a function of T. 



9 



FIG. 5. Configurational entropy Sconf- (a) Role of 
the different contributions in tfie estimate of Sconf- (b) 
i5/s-dependence of tfie configurational entropy for different 
p. Lines are extrapolation to lower Eis value based on the 
assumption of linear relation between Eis and l/T at low 
T. (c) Comparison between the Sconf values obtained using 
Eq.g and Eq.|l^ for p = 1.0 g/cw? . Symbols refer to different 
P{Eis{T),T) distributions, shifted to maximize the overlap 
for different T-values. 
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FIG. 6. Oxygen-oxygen pair correlation function, for equi- 
librated liquid configurations, and inherent structures. Upon 
lowering temperature, the first peak in both cases becomes 
sharper, and the intensity between the first and second peaks 
decreases. The smaller changes seen in the case of the inher- 
ent structures offers an estimate of that part of the change 
due to configurational change upon cooling, as opposed to 
thermal effects. 




FIG. 7. Pair correlation function for equilibrated liquid 
configurations, and inherent structures of the frequently stud- 
ied binary Lennard- Jones mixture between T = 0.44 and 
T = 5 in units of the pair potential energy depth [the model 
parameters are defined in Ref. Q]. Upon lowering tempera- 
ture, there is no apparent change in the structure of the IS, 
similar to the behavior we observe for water at large p. 
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FIG. 8. Average number of neighbors n(r) as a function of 
distance r. For all densities, an average of nearly four neigh- 
bors is reached at the first neighbor distance, showing the 
tendency for short range tetrahedral ordering, even at very 
high density. 
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FIG. 9. Distribution of the coordination number of 
molecules for both the instantaneous equilibrated configura- 
tions, and the inherent structures. In all cases, the fraction 
of four-coordinated molecules increases with decreasing T . 
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FIG. 10. Fraction of four coordinated water molecules: (a) 
equilibrated configurations and (b) inherent structures. 
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FIG. 11. Extrapolation for to lower T to estimate 

nu) = 1- 
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